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Abstract 

We consider supersymmetric model with Peccei-Quinn symmetry and study effects 
of saxion on the evolution of the universe, paying particular attention to the effects 
of thermal bath. The axion multiplet inevitably couples to colored particles, which 
induces various thermal effects. In particular, (i) saxion potential is deformed by 
thermal effects, and (ii) coherent oscillation of the saxion dissipates via the interaction 
with hot plasma. These may significantly affect the evolution of the saxion in the early 
universe. 



1 Introduction 



Peccei-Quinn (PQ) mechanism (TJ [2] provides an elegant solution to the strong CP problem. 
Promoting the ^-parameter to a dynamical variable, 9 = is realized after the spontaneous 
breaking of the PQ symmetry. In addition, in such a framework, a very light scalar field, 
called axion, shows up [3J H] . The coherent oscillation of the axion field is a viable candidate 
of dark matter of our universe if the PQ symmetry breaking scale is around 10 12 GeV [5]. 
Thus, the PQ mechanism has been attracted many attentions not only from particle-physics 
point of view but also from cosmology point of view. 

If we embed the PQ mechanism into the framework of supersymmetry (SUSY), which 
is a prominent candidate of the model beyond the standard model, serious cosmological 
difficulties may show up. In particular, because of the supersymmetry, there should exist 
superpartners of the axion, i.e., axino and saxion. Thermally produced axino and saxion, as 
well as coherent oscillation of saxion, may significantly affect the evolution of the universe. 
In order to construct cosmologically viable supersymmetric PQ model, it is important to 
understand how the axion multiplet affects the thermal history of the universe. 

The purpose of this paper is to study the effect of the saxion on cosmology, paying 
particular attention to the thermal effects on properties of the saxion. In hadronic axion 
models [SJ [7J, axion multiplet couples to extra colored multiplets, which include colored 
fermions (called PQ fermions), to solve the strong CP problem. Even in DFSZ-type model 
[HI [9] in which the ordinary Higgses have non- vanishing PQ charges, extra colored multiplets 
are likely to be introduced to avoid domain-wall problem. The saxion acquires thermal 
mass of the order of yT (with y being coupling constant) when the cosmic temperature T is 
higher than the mass of PQ fermions. Even if the temperature is lower so that the densities 
of the PQ fermions (and their scalar partners) are negligibly small, thermal effect induces 
logarithmic term in saxion potential [10J. These thermal effects shift the minimum of the 
potential, which may result in too large amplitude of the coherent oscillation of the saxion. 
Furthermore, the coherent oscillation of the saxion may dissipate via the interaction with hot 
QCD plasma. Taking these effects into account, we study cosmological history of a SUSY 
PQ model. 



2 Thermal Effects on Saxion 

In supersymmetric PQ model, there exists supermultiplet A (called "axion supermultiplet" ) 
which includes the axion field. Saxion (denoted as a) is a real scalar field in the axion 
multiplet: 

A = -=(a + ia) + V26a+ (F-term), (2.1) 
a/2 

where a and a are axion and axino, respectively. (Here and hereafter, the "hat" is for 
superfields.) The saxion potential is lifted only by the supersymmery breaking effects, and 
hence is very flat at T = 0. 
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The axion multiplet couples to colored multiplets, which include PQ fermions, to solve the 
strong CP problem. In this study, we consider the case that the masses of PQ fermions are 
related to the amplitude of a. The PQ fermions become massless when a takes a particular 
value. By integrating out PQ fermions, the axion multiplet couples to gauge multiplets (in 
particular, that of SU(3)c) as 

A„t = £|r J d 2 6AW a W a + h.c, (2.2) 

where VV is the field strength supermultiplet of SU(3)c- (Here and hereafter, a« = gf/4:ir, 
with gi being the gauge coupling constant of C/(l)y, SU(2)l, SU(3)c for i — 1, 2, and 3, 
respectively) 

The potential of the saxion is affected by particles in thermal bath (in particular, by 
those with SU(3)c quantum numbers). In addition, scattering processes may also become 
important. In the following, we discuss possible thermal effects on the saxion. 



2.1 Thermal Potential 

The first subject is thermal effects on the potential (more precisely, the free energy). The 
PQ fermions (and their superpartners), which we call Q, affect the potential as 

V T (a) = V {a)+g%W?\m Q {a))+gfv?\m Q {a)), (2.3) 

where Vq is the zero-temperature potential while Vj, and Vj, are contributions of bosonic 
and fermionic fields, respectively, and m,Q(a) is the mass of the fields which couple to o. 
(Extra contributions exist if particles other than PQ fermions couple to a.) At the one- loop 
level, the functions are given by [11] 



V 



rp4. poo 

(B/F) (M) = ±i- / dzz 2 In 1 T e-^M^] , (2.4) 

7T 2 Jq L J 



where upper and lower signs are for and V^ F \ respectively. In addition, g^^ are the 

number of degrees of freedom. (<7q = 1 for one complex scalar field, and g^ = 1 for one 
chiral fermion.) 

' and Vj. ' change their behavior at M ~ T. When M < T, Vj ' and Vj 1 are well 
approximated as 

V^\M) ~ -^T* + ^T 2 M 2 + 0(TM 3 ), (2.5) 

\4 F) (M) ~ _!^T 4 + lr 2 M 2 + 0(M 4 ), (2.6) 

and hence the so-called thermal mass shows up. When a satisfies the condition Wjq{(j) < T, 
Vr(cr) acquires a new term of ~ T 2 w?q{o\ For M > T, on the contrary, V^ and V^ 
rapidly go to zero. 
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If the temperature is high enough so that T 3> rriQ^ao), where o"o is the zero-temperature 
vacuum expectation value (VEV) of a, the PQ fermions are in thermal bath as relativistic 
particles even when a ~ a . In such a case, the expectation value of a is determined by the 
thermal mass term. Because Vj, {mQ{a)) and Vj. (rriQ^cr)) are minimized at a = oq (where 
oq satisfies mqipQ) = 0), we expect only one minimum of the potential at a ~ <jq at high 
enough temperature. The situation changes when the temperature is lower than m^fdo). In 
such a case, the number densities of Q are Boltzmann suppressed at a ~ a - Thermal effects 
are negligible around the zero-temperature minimum and there exist a minimum at a ~ a . 
(Effect of the logarithmic thermal correction to the potential [10J, which affects the position 
of the minimum, may be sizable. See the discussion below.) On the other hand, if Vo(cx) is 
flat enough, the minimum at a ~ <tq also remains. In such a case, there are two distinctive 
minima of the potential; even though there exists a minimum at o ~ a , the scalar field a 
can be thermally trapped in the different minimum close to the symmetry enhanced point. 
Then, the minimum at o ~ oq disappears when the temperature is so low that the shape of 
the potential around a ~ oq is determined by the zero-temperature potential. 

The saxion potential has another important term which is not taken into account in 
Eq. f)2.3p . Because the free energy of hot QCD plasma has a contribution proportional to 
gf (T)T 4 [12], and also because g%(T) depends on rngfa) (Le., the mass of colored particles) 
if rriQ^a) > T, the free energy has the following term [lO)^ 1 ! 

V L (a) = a L aj(T)T 4 \n\m Q (a)\ 2 , (2.7) 

where is a constant a bit larger than 1. 



2.2 Scattering Processes and Dissipation 

Next, we consider the scattering processes and dissipation. 

It has been known that scatterings of thermal particles contribute to the production 
processes of axion, axino and saxion. Define the yield variable as 



Y - nx 

Y X = : 



(2.* 



where nx is the number density of particle X while s = ^g*(T)T 3 is the entropy density 
(with g* being the effective number of massless degrees of freedom). Then, if the interaction 
of the axion multiplet with the minimal supersymmetric standard model (MSSM) particles 
is dominated by the operator given in Eq. (12. 2p . thermally produced axino abundance in 
radiation-dominated era is given by 



Y 



(th) 



mm 



RD 



r a (cq) ,0.20 x c^ln 



0.0977^ 



«3 



R 



10 7 GeV / V 10 11 GeV 



-2 



(2.9) 



# lr The logarithmic term is expected to disappear in the region where toq((t) is smaller than ~ T . In such 
a region, however, and dominate, so the effect of the logarithmic term is not important. 
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where Tr is the reheating temperature, and l^ cq ^ — 1.8 x 10 -3 is the thermal abundance of 
axino. A precise calculation of the abundance of thermally produced saxion is not available 
yet. Here, using the fact that the production processes of saxion and axino are governed 
by the same supersymmetric interaction, we approximate the thermally produced saxion 
abundance as [Yj th ^Rx> — |[^s ]rd- 

The thermally produced axino may be cosmologically harmful [H]. If it is the lightest 
superparticle (LSP) and is stable, it contributes to the present dark matter density. In such 
a case, a very stringent upper bound on the reheating temperature is obtained in order not to 
overproduce axino. (For recent discussion, see, for example, [T5| IT6] and references therein.) 
Even if the saxion is unstable, the LSP produced by the decay of axino contributes to the 
present mass density of the universe. If all the LSPs produced by the axino decay survive 
until today, low reheating temperature is needed to avoid the overproduction of the LSP. 
In the case of unstable axino, these cosmological difficulties can be avoided if the LSP has 
large enough pair annihilation cross section [17]. This may be the case if, for example, the 
LSP is neutral Wino [HI [191 [20] . in addition, i?-parity violation may be another possibility 
to avoid the cosmological difficulty. 

Existence of hot plasma not only increases but also reduces the number density of particles 
in the PQ sector. Such an effect is particularly important in studying the evolution of saxion 
oscillation. The dissipation rate of the saxion is related to the bulk viscosity of hot plasma 



HI [22], and is estimated aa^ 2 . 



9a§ T 3 



diss 



128tt 2 In c^ 1 F 2 : 



(2.10) 



where, if the saxion is far away from the minimum, F a should be replaced by the effective 
axion decay constant Fa° S ^ (which is ~ a). Notice that, if a interacts with particles other 
than MSSM gauge multiplets, T^ss receives contributions from those extra particles. We will 
see that this happens in some case. 

The dissipation rate Tdiss is an important quantity in studying the saxion oscillation in 
the early universe because Tdiss may become larger than the expansion rate of the universe 
H. For the case that Tdiss *C H, the saxion oscillates if the initial amplitude is non- vanishing. 
Then, with parabolic potential, for example, the energy density of the saxion decreases as a~ 3 
with the cosmic expansion (where a is the scale factor). On the contrary, if T^ss is sizable, 
we cannot neglect the effect of dissipation. In particular, if r diss > H, coherent oscillation 
of the saxion dissipates within a time scale shorter than the cosmic time. In radiation 
dominated universe with <7*(Tr) = 228.75, Tdiss given in Eq. (12.101) becomes larger than H 
whenT > 3xl0 5 GeV (5xl0 7 GeV, 7xl0 9 GeV, lxlO 12 GeV) for F a = 10 9 GeV (10 10 GeV, 
10 11 GeV, 10 12 GeV). One should note that, if r diss > H, yj th) becomes comparable to . 
Thus, in such a case, the particles in the axion multiplet are fully thermalized. 

Cosmology of models with saxion was considered in the framework in which the saxion 
field plays the role of the flaton field for thermal inflation [231 12^1 125] . We will consider a 



# 2 Thc authors thank K. Mukaida for helpful discussion on this issue. 
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different class of models, taking account of the thermal effects discussed above. We will see 
that the evolution of the saxion in the early universe can be significantly affected by these 
effects. 



3 Explicit Example 
3.1 Model 

To see the thermal effects on the evolution of the saxion, we consider the supersymmetric 
PQ model with the following superpotential: 

W = Wmssm + XS(XX - f) + yXQ Q + , (3.1) 

where Wmssm is the superpotential of the MSSM, and X (+1), X (—1), Qo (0), Q+ (+1), 
and S (0) are chiral superfields. (We denote the [/(l)pQ-charge in the parentheses.) Qo 
and Q + are in the 3 and 3 representations of SU(3)c', here one pair of PQ fermions are 
introduced. One combination of X and X becomes the axion sup er mult iplet. In addition, 
y, A, and / are constants which are chosen to be real and positive. 

In the minimum of the supersymmetric scalar potential, we obtain the relation XX = f 2 . 
Such a constraint eliminates one combination of X and X from the low-energy spectrum 
and X and X are decomposed into the axion multiplet A and heavy multiplet which we call 
X . Without SUSY breaking terms, the relative size of X and X is undetermined; such a 



flat direction corresponds to A. The SUSY breaking terms relevant for our study ara^ 3 

V soit = m\\X\ 2 + m\\X\ 2 + m 2 s \S\ 2 , (3.2) 

where we expect that m 1; m 2 , and m$ are of the order of the gravitino mass m 3 / 2 . With the 
SUSY breaking terms, the VEVs of scalar fields are fixed. 

At the minimum of the potential, the axion multiplet A and the heavy multiplet X are 
embedded into X and X as 

X = F a cos f3x — A cos f3x + X sin f3x, (3.3) 

X = F a sin p x + A sin f3 x + X cos (5 X , (3.4) 

where F a is the axion decay constant which is given by 

*-<*>• + <*>' = 5^. (3-5) 

with 

2 2 

mi — mi 

cos2/3 x =e = ^ f- 3.6 

m\ + m 2 



^ 3 The tri- linear interaction SXX and linear term in S may also exist in V so ft, with which the axino mass 
is generated. In the present model, the axino mass is at most -j=m a . 
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By integrating out Qo and Q + , we obtain the interaction given in Eq. ( 12. 2p . Then, the saxion 
mass is m a = 2mim 2 / \Jm\ + m§. 

We note here that, in the present setup, the decay rate of saxion r CT is governed by the 
£ parameter defined in Eq. (I3.6p . In hadronic axion model in which axion multiplet does 
not couple to the Higgs multiplet, a dominantly decays into axion pair; decay rate of such 
a process is given by 

e < (3.7) 



64vr F 2 



The saxion also decays into gauge-boson pairs (in particular, the gluon pair) via the interac- 
tion given in Eq. (12. 2p . However, the saxion-gluon-gluon interaction is loop suppressed and 
hence such a process is expected to be subdominant (as far as £ ~ 1); indeed, the decay rate 
for the process cr — >■ gg (with g being the gluon) is given by 



V ^"~ 2^fV (3 ' 8) 



a 



Numerically, T a _, aa > T a ^ gg when f > 0.05. 

X or X may directly couple to the MSSM Higgses if they have non-vanishing PQ charges, 
which gives the decay mode of saxion into the Higgs boson pair. For example, we may 
introduce the following term into the superpotential: 



W X ^H d = j^-X 2 H u H dl (3.9) 

where H u and H d are up- and down-type Higgses, respectively. Using the fact that the 
so-called /i-parameter is generated by the VEV of X, we obtain the decay rate as 

i 3 / 2 \ 4 / a\ 1/2 

1 m a ( n \ ( m h x 



r -=^l^J ■ (310) 

Br (a — > hh) depends on various parameters. If fi ~ 0(m a ), Br(a — > hh) and Br(a — > aa) 
are of the same order. In models with smaller m a (like the gauge mediation), the decay 
process of a — > hh is kinematically blocked. 

Since we are interested in the case of F a ~ 10 9-12 GeV, the lifetime of the saxion becomes 
very long. If there exists relic saxion in the early universe, the decay of saxion produces axion 
which survives until today and behaves as an extra radiation component (so-called "dark 
radiation" ) [23J . We will study implications of such a relativistic axion in the present model. 



# 4 If X, instead of X, couples to the Higgses, the potential at X ^> f may be modified due to extra thermal 
effect, which may change the following discussion. We do not consider such a case. 
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3.2 Thermal Potential 



Now, we are at the position to discuss the potential of the saxion. We first comment on 
the so-called Hubble-induced mass. In the inflaton dominated era before reheating, Planck 
suppressed interactions of X and X may induce their effective masses of the order of the 
expansion rate of the universe. Even in radiation dominated universe, this may be the case 
[26] . (We call such masses as Hubble-induced masses.) The sizes of the Hubble-induced 



masses crucially depend on the model, and we simply parametrize the potential aa# 5 



V T = m 2 |X| 2 + ^ + 2A^vfV7W) + ^^ 

+2l4 B) (AvW+/W) + 2l4 F) (AvW+7W)- (3-11) 

Here, fhi are effective masses which may include the effect of Hubble induced masses; fhi ~ rrii 
when H <C 7713/2, while fhi ~ H if the Hubble induced masses exist and dominate. In Eq. 
(13. lip . N c = 3 is the color factor. In addition, in our numerical study, we take — 1. 

Once the amplitudes of X and X become non-vanishing, Q and Q + become massive 
with the mass of y\X\. In addition, 5* and X acquire the mass of Aa/|X| 2 + |X| 2 . After 
the PQ symmetry breaking, the product XX is fixed to be f 2 + 0(yTfhi). Thus in the 
situation of our interest, we can eliminate X (or X) from the potential. We choose X as the 
independent variable because X plays the role of the axion multiplet when X ^> X, which 
is the case in the following discussion. 

Although S and X also affect the thermal potential, it is instructive to consider the case 
that the thermal effects are dominated by those of Qo and Q + . In such a case, the potential 
is approximated as 

Vt ~ rh\\X\ 2 + ^ + a T y 2 T 2 j^0(T - yf/\X\) - a L a\T' In |A| 2 , (3.12) 

where ax is a constant of 0(1). When the temperature is high enough, the potential have a 
minimum at X ~ xt, where 




XT = J^-f. (3.13) 



(Hereafter, we call the minimum at X ~ xt as "trapping minimum.") The trapping min- 
imum disappears when T becomes smaller than the effective mass of the PQ fermions at 
X ~ xt- This happens when the cosmic temperature is ~ T c , where T c is the solution of the 
following equation: 

T c = ^j— (3.14) 



# 5 If X 01 X couples to extra fields, contributions of those extra fields should be taken into account. This 
is the case if, for example, Qo and are embedded into complete multiplets of the gauge group of grand 
unification. 
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(Notice that, when the Hubble-induced mass is sizable, fh\ depends on background temper- 
ature.) On the contrary, at low temperature, the T 2 term in the potential does not exist 
at the region with X <C yf 2 /T. If so, there may exist a minimum at X ~ xl, where xl is 
given by 



2 _ a L a 2 T* + V(a L alT*y + jgfgfg 
x L = ^ . (3.15) 



When aLu\T A > mim 2 / 2 , xl ~ a 1 / 2 'a^T 2 jm\, which can be significantly larger than /. 
Behavior of the potential at X ~ xl depends on the effective mass of the PQ fermion at 
such a field value. If the PQ fermion mass is smaller than T, Q and Q + (and their scalar 
partners) are thermalized as relativistic particles at X ~ xl, and hence the T 2 term in the 
potential exists. In such a case, X ~ xl is not the minimum of Vr, but there exists only one 
minimum at X ~ xt- This is the case when T > T*, where is obtained by solving the 
following equation: 

vf 2 

T* = JLI. (3.16) 

(In defining T*, effects of S and X are not included; their effects will be separately considered 
later). If the Hubble-induced mass is negligible, T* = a 3 1 ^ 3 y 1 ^ 3 m\^ 3 f 2 ^ 3 (where we have 
assumed a 3 ?/ 2 (/ /m 3 / 2 ) 5 > 1, which is the case in the situation of our interest.) On the 
contrary, for T < T*, the PQ fermions are heavier than T at X ~ xl, and we expect that 
the potential has a minimum at X ~ x^. 

In summary, the thermal potential behaves as follows: 

• T < T c : There is only one minimum at X ~ x^. 

• T c < T < T^: There are two minima at X ~ xl and X ~ xt- 

• T > T*: There is only one minimum at X ~ xt- 

We should also comment here that, if xl(T c ) ~ 0(f ), two minima become indistinguishable 
at T ~ T c . In such a case, trapping minimum, which exists at high temperature, smoothly 
merges to the zero temperature minimum. Our numerical calculation indicates that this 
happens when y < lyrn^JJ ■ In this case, the coherent oscillation of the saxion is not 
induced. 

In Fig. CD, we plot Vr given in Eq. (13. lip for the temperature around T ~ T c , neglecting 
the effects of S and X. (Here, we take fhi = rrii and the effects of Hubble-induced masses 
are omitted.) As one can see, when the temperature is 0(T C ), the trapping minimum disap- 
pears. In addition, two distinct minima exist just before the disappearance of the trapping 
minimum. Behaviors shown in the figure are consistent with the previous discussion. 

So far, we have considered the case that the thermal effects from S and X are negligible. 
However, if Xf < T c , the potential around X ~ / is deformed. In particular, if A is sizable, 
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X / f X / f 

Figure 1: Thermal potential as a function of Xj f. We take / = 10 10 GeV, mi = m,2 = 1 TeV, 
and y = l. The solid (dotted, dashed) line is for T = 0.69T C , (T = 0.67T C , T = 0.65T C ). We added 
a constant to the potential. In addition, the vertical axis is normalized by the zero-temperature 
expectation value of the potential (Vo) = m\(X) 2 + m^X) 2 . The left figure is for the case that the 
contributions of S and X are negligible, while the right figure is for the case with A = 0.01. 

there exists a minimum at X ~ / at T ~ T*. In Fig. [TJ we also show one of the examples 
of such a case, taking A = O.Oly. Such an extra minimum may affect the evolution of the 
saxion, in particular the trapping process. (See the following discussion.) 

The thermal trapping of the saxion may be cosmologically important, although its im- 
plication strongly depends on the details of the model. In the following, we discuss possible 
effects of the thermal trapping in the present model. We pay particular attention to the en- 
ergy density of axion produced from the saxion oscillation. If the saxion field is once trapped 
in the trapping minimum, it starts to oscillate with the amplitude of ~ xt at T ~ T c . Then, 
such a saxion oscillation decays with long lifetime. If a dominantly decays into axion pair, 
the produced axion survives until today and behaves as dark radiation whose abundance is 
constrained from the cosmic microwave background (CMB) anisotropy and big-bang nucle- 
osynthesis (BBN). 

3.3 Case with y < A 

Thermal history depends on the mass spectrum of particles in the PQ sector. We first 
consider the case with y < A; in such a case, S and X are always heavier than Qq and Q + . 
Then, at T < A/, we neglect the thermal effects from 5" and X. The cosmological implication 
of the saxion depends also on its interactions. Thus, for simplicity, we mostly concentrate on 
the case that the interaction of the saxion with the Higgses given in Eq. (13.91) is negligible. 

We start our discussion with studying the behavior of T*, which depends on the size of the 
Hubble-induced mass. In the following study, we consider the case that the Hubble-induced 
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10 7 GeV, y = 0.01 



10 7 GeV, y = 1 



10 12 - 



T 




10 3 
m a (GeV) 



10 4 10 2 



10 3 
m ff (GeV) 



10 4 



Figure 2: The dotted lines are contours of constant T*. The solid lines are contours of constant 
R a / r (see Eq. (|3.2ip ) for £ ~ 1. (Numbers in the figures are the values of R a / r .) In the darkly- 
shaded (dark blue) region, AN V becomes smaller than 1 if Fa/T a = 1CP 6 . In the lightly-shaded 
(light blue) region, AN U becomes smaller than 1 if the axino decays just before the BBN. Above 
the dot-dashed line, Eq. (13.16P does not have solution. Below the dot-dot-dashed line, Tdi ss > H is 
realized at a certain epoch after T = T C . Here, Tr = 10 7 GeV, and the Yukawa coupling constant 
is taken to be y = 0.01 (left) and 1 (right). 



mass is sizable and approximate 

fhi = max(mj, H). (3-17) 

The Hubble-induced mass relaxes the constraint compared to the case with fhi = rrii. 

In Figs. [2] and [3], we plot the contours of constant T* on m a vs. F a plane for fixed values of 
y and Tr. (Here and hereafter, we approximate / ~ F a and m\ ~ m 2 ~ fn a .) The behavior 
of T* can be understood as follows. For small enough F a , effect of the Hubble-induced mass 
is irrelevant in solving Eq. ( 13.161) . and T* is independent of Tr. On the contrary, for larger 
value of F a , effect of the Hubble-induced mass becomes important. Let us consider how T* 
behaves when fhi ~ H . In the radiation-dominated universe, H is proportional to T 2 , and 
hence xl — a 1 / 2 a^T 2 / H becomes insensitive to T. In such a case, T*, which is the solution of 
Eq. (I3.16p . becomes independent of mi. On the contrary, in the inflaton dominated universe, 
H oc a~ 3//2 while T oc a -3 / 8 , and hence 

H(T > Ik) ~ H(T R ), (3.18) 

where T is the temperature of dilute plasma produced by the decay of inflaton. Then, xf} 
is proportional to T 2 as the temperature of dilute plasma increases, and Eq. (13.161) does 
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not have a solution if F a is too large. Consequently, the thermal trapping of the saxion 
is not guaranteed if F a > gl 3 ^ 16 y^ 1 ^ 2 a^ 2 m!/ 8 T^ 4: Mp l 8 . We conservatively assume that no 
constraint is obtained in such a case. 

If the saxion is once trapped in the trapping minimum, the saxion starts to oscillate with 
the amplitude of 0(xt) when the trapping minimum disappears. Even if X and X are of 
0(f) just after the PQ phase transition, we expect that the thermal trapping occurs if T* 
is lower than the maximum temperature of the universe T max ~ H^Mp^T^ 2 , with H ini 
being the expansion rate during inflation |27j . 

In Figs. [2] and El we also show the region in which the dissipation rate becomes larger than 
H at a certain cosmic temperature below T c , where the dissipation rate is evaluated with 
Eq. ( 12.1 Op by replacing F a with the amplitude of X. In such a region, even if the trapping 
happens, the saxion oscillation dissipates away. Then, the relic saxion is dominated by 
thermally produced one. Notice that the region with y < l-sjmij ' j ', in which two minima 
merge smoothly, is mostly covered by the region with T^iss > H. One may think that the 
scattering processes among the particles in the PQ sector also contribute to the dissipation 
process. At T < T c , only the particles in the axion multiplet A (i.e., axion, axino, and 
saxion) can contribute to the dissipation process because the masses of other particles in the 
PQ sector are larger than T c . The axion multiplet particles are thermalized only when Tdiss 
becomes comparable to H. In addition, the interaction of the axion multiplet is suppressed 
by inverse powers of Fjf s ^ because the axion is a Nambu-Goldstone boson. Dissipation rate 
due to the scattering processes among the PQ sector particles is estimated to be ~ j- ^ )4 

and is smaller than T d „, given in E q . « (with F a -» F<*>) in the parameter reg'ion of 
our interest. 
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If the dissipation of the saxion oscillation is negligible, the saxion oscillation with large 
initial amplitude may result in the overproduction of the axion dark radiation. To see when 
there may exist such a problem, we calculate the density of the oscillating saxion. The yield 
variable of the saxion is estimated as 

yjosc) „ 45 m a| 



7T 



where <7r is the saxion amplitude at T = Tr, and m a is the effective mass of saxion which may 
include the effect of the Hubble-induced mass. In our numerical analysis, we use g*(T R ) = 
228.75. (If the saxion starts to oscillate after the reheating, T R should be replaced by 
the temperature at which the oscillation starts, and ctr should be identified as the initial 
amplitude of the saxion oscillation.) If the oscillation of the saxion starts before the reheating, 
^o-(2r) is related to the initial number density using the fact that a oc T~ 3 / 8 during the 
inflaton-dominated era. Denoting the temperature at which the saxion starts to oscillate as 
Ti (with Ti > T R ), 

n a (T R ) ~ (^) n a (Ti). (3.20) 



T 

To see how large the saxion abundance is, we calculate the following ratio 

E> Po-(^dcc) /q qi\ 

K a/r = 77 y 

Pr \J^dec ) 

where p a (tdec) and p r (tdec) are energy densities of the saxion and the radiation from the 
inflaton decay at the time just before the saxion decay, respectively. If R a / r > 1, the saxion 
dominates the universe. In Figs. [2] and [31 we show contours of constant R a / r . We can see 
that the saxion dominance occurs in large fraction of the parameter space. Then, when it 
decays, some amount of energetic axion is produced. 

The present density of the axion from the saxion decay is parametrized by using the 
effective number of extra neutrinos as 

AN V = N<& - 3 ee 3pa(tnow) . (3.22) 

Here, p a (t now ) and p u (t now ) are energy densities of the axion and neutrinos in the present 
universe, respectively, and the factor of 3 in the numerator is the number of generations 
of neutrinos. From the latest WMAP observation (28], the best-fit value of the effective 
number of neutrinos (which includes the standard-model contribution) is obtained as Njf S ^ = 
4.34l°;| (68% C.L.). In addition, the analysis of 4 He mass fraction generated by the BBN 

reactions indicates n!? s) = 3.68±g;fg (2a) |2S]. 

If the saxion abundance is too large, AN U exceeds ~ 1, which conflicts with observations. 
In estimating AN,,, we should notice that the entropy production due to the axino decay 
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may be sizable in the present case. AN V becomes smaller as the decay rate of axino Ta 
gets smaller. The axino decays into gaugino and gauge boson pair and the total decay rate 
of a is smaller than that of o if £ is close to 1. Typically, Ta/T a ~ 0(-/Vjof /2tt 2 £ 2 ), where 
Ni is the number of final states and i depends on the dominant decay mode; if the axino 
mass is well above the gluino mass, Ni = 8 and on = 013. However, if the masses of the 
axino and the gauginos are degenerate, T s becomes suppressed. So the ratio Ta/T a is model 
dependent. Notice that, if T s is too small, the axino decays after the BBN, which is likely 
to spoil the success of the standard BBN scenario. We vary T s and estimate AN U . (We use 
the instantaneous decay approximation.) Here, the axino abundance is evaluated using Eq. 
f 12 .91) . and the axino mass is taken to be \m a . In Figs. [2] and El we show the region in which 
AN V becomes smaller than 1 if Ta/T a = 10 -6 , and also the region in which AN V becomes 
smaller than 1 if the axino decays just before the BBN. If Ta/T a = 10~ 6 , AN V > 1 for 
R a / r > 1 even with the entropy production due to the axino decay. If the axino decays just 
before the BBN, the effect of the entropy production is more significant. Here, the energy 
density stored in the axino sector is assumed to be fully converted to that of radiation after 
the axino decay. This may be due to the pair annihilation of the LSP or due to the decay of 
the LSP via i?-parity violation. Even with such an extreme assumption, AN V may be larger 
than ~ 1 if Br [a — y aa) ~ 1 in significant fraction of the parameter space. 

Such a problem may be avoided if the SUSY breaking parameters are tuned so that the 
£ parameter becomes relatively small. Indeed, if £ < 0.05, Br (a —y gg) becomes comparable 
to or larger than Br (a —y aa), and the production of the axion is suppressed. In the region 
with Ra/r > 1 of the figure, saxion once dominates the universe. Even so, AN V < 1 is 
possible if Br(o — y aa) ~ O(0. 1). In fact, as one can see, the constraints on Nu from 
the CMB anisotropy and BBN indicate non- vanishing value of AN U ; AN U ~ 1 is preferred. 
If Br(a — y aa) ~ 0.25, AN V = 1 is realized (with g* = 100 at the time of saxion decay). 
Another solution is to introduce the superpotential interaction given in Eq. (I3.9p . With 
such an interaction, a may dominantly decay into Higgs boson pair, and Br (a — y aa) can 



be O(0.1) j£3 In such a case, the saxion may also decay into the Higgsino pair, which may 
result in the overproduction of relic LSP. However, such a difficulty can be avoided with 
large enough pair annihilation cross section of the LSP. This is the case, for example, if the 
LSP is the neutral Wino [18]. If the saxion dominantly decays into MSSM particles, one 
should take account of the entropy production due to the saxion decay. The dilution factor 

cr/r 

3.4 Case with y > A 

Next, we consider the case with y > A. In such a case, if X ~ X ~ /, Qq and Q + become 
heavier than S and X. Even so, the trapping may happen if T* < A/; in such a case, 
discussion in the previous section holds. 



# 6 The interaction with the Higgses may enhance the dissipation rate, which may also help to reduce the 
density of the axion. The study of the dissipation rate in such a case will be given elsewhere [3U1 • 
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If T* > A/, the situation is rather complicated. This is because S and X may be fully 
thermalized as relativistic particles at T < T c . Then, the scattering processes with S and 
X significantly contributes to the dissipation of the oscillation of the saxion. In addition, if 
T* > A/, there may exist a minimum at X ~ / caused by the thermal effects of S and X, 
with which the trapping process to the trapping minimum may not happen. (See Fig. [TJ) 
The S and X can be produced from the scattering of particles in thermal bath (like gluon). 
Even though the interaction between the MSSM sector and the PQ sector (which consists 

of X, X and S) may be too weak to equate the temperatures of two sectors, we expect that 
the interactions among X, X and S are strong enough to thermalize the PQ sector; such a 
thermalization occurs through the superpotential interaction of XSXX. Thus, if the masses 
of S and X are smaller than T* (at X ~ /), constraint obtained in the previous section may 
not be applicable. This is the case when A < y 1 ' 3 (m a -/F a ) 1 ' 3 . 

3.5 Saxion from Scattering of Thermal Particles 

Finally, we comment on the effects of thermally produced saxion. As we have mentioned, 
scattering processes of thermal particles produce saxion. In particular, if T diss ~ H is realized 
at some epoch, the saxion are thermalized. In such a case, the overproduction of the axion 
may occur if the saxion dominantly decays into axion pair. However, the abundance of 
thermally produced axino is comparable to that of saxion. Then, entropy production due to 
the decay of axino (as well as due to the subsequent pair annihilation of the LSP) dilutes 
the axion abundance. If the decay rate of the axino is smaller than that of saxion, the effect 
of the dilution is large enough to make AN U smaller than ~ 1. 

4 Discussion 

In this paper, we have studied effects of saxion on the evolution of the universe paying 
particular attention to thermal effects on the saxion. Because the axion multiplet necessarily 
couples to colored particles (i.e., PQ fermions), saxion potential is inevitably deformed at 
high temperature. In addition, the coherent oscillation of the saxion may dissipate via the 
interaction with hot plasma. 

Taking account of these effects, we have studied the constraint on the supersymmetric 
PQ model with the superpotential given in Eq. (13. ip . Because of the thermal effect on 
the potential, the saxion field may be trapped in a false minimum before the start of the 
oscillation, which may result in an overproduction of relativistic axion. The constraint 
strongly depends on the model parameters. We have seen that, if no PQ sector particle 
remains in thermal bath after the onset of the oscillation, AN U ^> 1 may happen if Br (a — > 
aa) ~ 1. In addition, if the dissipation rate is large enough, the saxion oscillation dissipates 
soon after the start of the oscillation. 

In this study, we have concentrated on the model with the superpotential given in Eq. 
( 13. ip . However, the thermal effects may have significant effects in other classes of models. 
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